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Abstract 

Many popular models of fractionalized spin liquids contain neutral fermionic spinon excitations 
on a Fermi surface, carrying unit charges under a compact U(l) gauge force. We argue that 
instanton effects generically render such states unstable to confinement in two spatial dimensions, 
so that all elementary excitations are gauge neutral, and there is no spinon Fermi surface. Similar 
results are expected to apply to SU(2) spin liquids. However, fractionalized states can appear when 
the gauge symmetry is broken down to a discrete subgroup by the Higgs mechanism. Our argument 
generalizes earlier results on confinement in the pure gauge theory, and on the instability of the 
U(l) 'staggered flux' and 'algebraic' spin liquids with a Dirac spectrum for fermionic spinons. 
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I. INTRODUCTION 
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The anomalous normal state of the cuprate superconductors has focused a great deal of 

ization. One model, which has e njoy ed a 
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attention on exotic states with electron fractiona 

QSQQBBQ 



deal of popularity 

21 



22j,|23(, has the electron spin reside on neutral S = 1/2 fermionic 'spinons' which form a 
Fermi surface. The microscopic theory of the formation of spinons shows that each spinon 
carries unit charge of a compact gauge group, which is usually 
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2l\ U(l); however, somewhat more involved theories with 
a SU(2) gauge group have also been considered [3|, [16, 000. The spinons interact 
with each other, and also with other possible 'holon' excitations, via exchange of the gauge 
boson. We will focus here on the U(l) case, but our arguments appear to have a direct 
generalization to the SU(2) case. 



Polyakov 



24j argued many years ago that the pure compact U(l) gauge theory is al- 



ways confining at zero temperature in two spatial dimensions. In other words, there is no 
'Coulomb' phase with a massless 'photon' excitation, and oppositely charged static test par- 
ticles experience a confining potential which increases linearly with the distance between 
them. The confinement was induced by the proliferation of 'instanton' tunneling events. 
Given this fundamental result, one might wonder if similar effects are important in the 
models of electron fractionalization Q]. Indeed, in theories of quantum antiferromagnets 
with collinear spin correlations, such effects were crucial in disrupting possible spin liquid 
states and leading to ground states with bond order and confined spinons 0|. 

Many models of spin liquids have gapless fermionic spinon excitations, and in these cases 
the argument for confinement by proliferation of instantons is more delicate. The 'stag- 
gered flux' and related 'algebraic spin liquids' possess spinons with a relativistic Dirac-like 
spectrum, and consequently the instantons interact with an action which depends logarith- 



mically on their spacetime separation |2 
of a single instanton, several authors |lj 



B y ex amining the renormalization of the fugacity 
l3- have argued that this logarithmic coupling 
leads to a regime in which instantons are suppressed, and so a fractionalized spin liquid 
regime is stable. Similar arguments have also been advanced recently for cases where there 



are relativistic massless Bose fields 



28|. However, these arguments neglect the screening of 



the instanton-instanton coupling by bound instanton pairs at shorter scales. A number of 
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other authors [29|, |3Qj, |3jJ, |32| have mentioned this screening, and argued that it disrupts 
the binding of instantons in 2+1 spacetime dimensions; the instantons are always ultimately 
in a plasma phase, interacting via an exponentially screened interaction, and the spinons 
carrying dual 'electric' charges are consequently permanently confined (these arguments are 
presented here in Appendix |XJ) . Recently, the case for permanent confinement was put on a 
firmer footing by two of us [33]: an explicit duality mapping of a compact U(l) gauge theory 
minimally coupled to relativistic fermions showed that the plasma state of instantons was 
the only stable phase, and so there were no 'electrically' charged fermionic spinons. This 
result implies that the staggered flux and algebraic spin liquids do not exist in two spatial 
dimensions. 

The problem addressed in this paper is the case of a compact U(l) gauge theory coupled to 



non-relativistic matter, with the fermionic spinOT^fOTmingja ivnm sur 
this case has 
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ace. As noted above 



12 



13 



14. 



15 



as 



}een the focus of much attention 
20l l2l| in models of the normal state of the cuprates. For this case, Nagaosa [10] has 
advanced arguments that the dissipation caused by the spinon Fermi surface can suppress 



the instanton quantum fluctuations, and so stabilize the spin liquid state. However, using 



the method of Ref. 
between instantons 



33l we show below that Nagaosa's argument also neglects the screening 



34j . and that the spinon Fermi surface is generically unstable in two 



spatial dimensions at zero temperature. 

A surviving route to fractionalization in two spatial dimensions arises in physical situ- 
ations with a Higgs field which does not belong to the fundamental representation of the 




In the spin liquid case, such a Higgs field is provided by the pairing of 



gauge group 
the spinons 



| in a BCS-like state, and the formation of such a BCS condensate 
breaks the U(l) gauge symmetry down to Z 2 . The resulting Z 2 gauge theory does permit 
a deconfined phase, with spinons which carry only Z 2 gauge charges. Note, however, that 
because of the pairing, there is no spinon Fermi surface in this case either. 



42|, which allows 



For completeness, we also mention the 'chiral spin liquid' 0, Q, 
fractionalization in the presence of time-reversal symmetry breaking, but again without a 
spinon Fermi surface. 

This paper will only be concerned with matter in the fundamental representation of the 
U(l) gauge group. Our analysis below shows that the only possible exceptions to conven- 
tional confinement in such situations arise at special critical points where a certain relevant 
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coupling in the action accidentally vanishes: at least one parameter has to be tuned to a 
particular value for this to happen. At such critical points, the gauge theory is not in a con- 
ventional deconfmed Coulomb phase either, and it is likely that there are no quasiparticle 
excitations. These points are analogous to the critical Rokhsar-Kivelson points in quantum 



dimer models 
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Our analysis of the compact U(l) gauge theory coupled to fermionic spinons on a Fermi 
surface begins in Section |H] with a duality mapping to a sine-Gordon-like theory. The renor- 
malization group analysis of the latter theory is presented in Section ITTT1 and Appendix [Bj 
Appendix |X] contains a discussion of similar physics expressed directly in terms of the sta- 
tistical mechanics of instantons. 



II. DUALITY MAPPING 

A fundamental characteristic of the proposed spin liquid state j4j, is the overdamped 
nature of the transverse gauge field propagator. This damping arises from the low energy 
fluctuations of the spinon Fermi surface. We denote the components of the compact U(l) 
gauge field by a M with \i = 0,x,y, the zero'th component representing the imaginary time 
direction; also by dj we denote only the spatial components of a M . Then, after integrating 
out the fermionic spinon fields, the action, S, for^a^ at small wavevectors (k) and small 
frequencies (u) has the following singular form 



S 



d 2 kduo 
8tt 3 
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Here cr(k) is the spinon "conductivity", \d is the spinon diamagnetic susceptibility, \ is 
the spinon compressibility, and 7 is a damping coefficient characteristic of the spinon Fermi 
surface. All of the terms in S are characteristic of a fermionic spinon system with a finite 
density of states at the Fermi level (for the relativistic case with a Dirac spinon spectrum, 



there is a vanishing density o: 



field action quite different j33|). For a(k), it is conventional to assume |4 

, £ for H < 1 
a(k) ~ (2) 
1 /k for kl > 1 

where i is the mean-free path associated with scattering off static impurities. 



states at the Fermi level, and this makes the effective gauge 
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Our subsequent analysis is aided by writing S in a form which is explicitly gauge invariant. 
We define the electromagnetic field F^ u = p^a u — Pva^ where po = u, Pi = fcj. Then, it is 
easy to see that 



d 2 kdu 



FioFjo 



(3) 



4 V k 2 e 2 J 13 2 [ e 2 W \ k 

The expression © for S also includes a regular F 2 t/ /(4e 2 ) term in the action for the gauge 
field; such a term will invariably be generated by integrating out high energy fluctuations of 
matter fields at short scales. 

We may now proceed to dualize the above action for the gauge fields. Before doing so, 
it is necessary to explicitly account for the compact nature of the U(l) gauge group on 
the scale of underlying lattice. We do this by discretizing space and time, and writing the 
compactified lattice version of the action (jHJ) as 
1 



s= E 



X,y,T,T' 



-(Fij(x, t) - 2im ij (x, T))V t {x - x , r - T l )(F ij (x l , t') - 27m ij ( y x', r')) 



(4) 



+-{Ai{F i0 {x, r) - 2nn i0 (x, T)))V t {x - x' , r - r')(Aj(Fj (x' , r') - 27m j0 (x', r'))) 



^ e X,T 

where the transverse and the longitudinal parts of the interaction in the Fourier space read 

(5) 



V t (k,u) = \y +Xd + -~: 

e 2 



k 2 



(6) 



The integers n^ u (x, r) serve to account for the compact nature of the gauge fields, in the spirit 



of the Villain approximation. On a lattice, the electromagnetic tensor F. 



fiv '-^iva^-apcj'—^p 



A p a a , 



where is the lattice (discrete) derivative. 

Performing the Hubbard-Stratonovich transformation using fields c^ u residing on the 
plaquettes of the lattice, we find 



s= E 



x, x' ,t,t' 



1 



Cij(x, r)V t (x -x',t - T')cij(x', t') + -e c i0 (x, t)< 5 i:j 5 xx '5. 



+ 



|A| 2 Ve 2 



AiAAc^x'y) 



5xx>$tt' + Vi(x -x',t -t')\A\ 2 

+^E c ^( x ' r )( i7 A^( x ' r ) ~ 27rn^(x,r)) 



(7) 



The integral over the gauge fields can now be performed exactly, and it enforces the con- 
straint 

Cp,(x,T) = -e Mm A a $(x,r). (8) 

The sum over the integers n M! ^(x,r) can also be performed, and it forces the new variable 
<&(x,t) to be an integer as well. From the constraint in Eq. (JBJ) it immediately follows, 



A,;C,;nfx,r) 



0. 



(9) 



and the corresponding part in Eq. 0, which includes the longitudinal part of the interaction 
Vi, completely drops out. Relaxing the integer constraint on $ by introducing a small 
fugacity y in the usual way 4^|, upon return to the continuum limit we find the dual sine- 
Gordon action for the instantons 



S, 



G 



d 2 kduo 



y\uj\a(k)/k 2 + Xd + l/e 2 
2y J d 2 xdr cos(27r$(x, t)). 



+ e z k 



2/„2 



(10) 



28. The 



A corresponding sine-Gordon theory for the relativistic matter case appears in Ref. 
problem of confinement in the original theory now reduces to understanding the possible 
phases of S s q, to which we turn in the next section 



III. RENORMALIZATION GROUP ANALYSIS 

To be specific, we assume that in the sine-Gordon theory the wavevectors \k\ < A, 
with A£ < 1, and frequencies |u;| < Q. To simplify the calculation, we also take Q 3> 
(Xd + e~ 2 )A 2 /£, so one can write the sine-Gordon field propagator for low momenta and 
frequencies as 

G-\k,u) = + a k k 2 + a^ 2 , (11) 

where the bare (unrenormalized) coefficients are a k = e 2 , and = 0. Our methods below 
are easily generalized to the clean limit case where cr(k) ~ l//c, and the results are very 
similar. 

Let us integrate out the short-distance modes with A/b < \k\ < A, and < Q, with 
ln(6) 1. The parameters for the remaining long-distance modes are then changed as 

m = i, (12) 
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a k (b) = a k + I 2/ 2 e -G > (^=o,r=o) J ^g dT \^ e G>(s,r) _ 1); (13) 

a w (6) = b 2 [a u + y 2 e -G>(x=0,r=0) | ^ d ^ e G>(S,r) _ ^ ( M) 



y (6) = 6 2 ye-i G >^ =0 ' T= °), (15) 
where the correlation function is defined as 

G>(x,t) = f p.e***"G$,u,). (16) 

J\/b<\k\<A {2tt) z J-n 2tt 

A crucial feature of the above recursion relations is that although initially the coefficient 
= 0, for b > 1 a w (b) > 0. This expresses a simple physical effect that the interaction 
between two distant charges in three spacetime dimensions in a medium with a finite polar- 



izability is always Coulombic (~ 1/y/x 2 + r 2 ) in three spacetime dimensions The only 
apparent way the generation of this term can be avoided would be if the propagator G(k, uS) 
would be completely independent of frequency, in which case G > (x, r) ~ ln(6)e _ ' z ' A 5(r), ren- 
dering the integral in a w (b) to be exactly zero. This unphysical situation would correspond 
to an essentially two dimensional sine-Gordon theory, which would then have the standard 
Kosterlitz-Thouless transition. 

In the present case, on the other hand, we find 



dd k y 



2 



d\n(b) 2nd k 



+o(n (17) 



2^ + ^ + 0(f), (18) 



dln(6) 4a k 

where we have introduced the dimensionless combinations y = y/A 2 Q, d k = a k /Q, and 
duj = a^Q/A 2 , and introduced a smooth cutoff as described in Appendix iBl 

As announced earlier, a finite coefficient becomes generated to the second order in 
fugacity, and then becomes a relevant coupling. With this term included the fugacity flows 
according to 

T^7T = (2 ; 1 )y + Q{f), (19) 

dHb) Svr^ - (l/2£)^ 

and always becomes relevant at long length scales. In Eq. ([T9*)l we assumed dkd^ ^> (l/2£) 2 ; 
other limits and more complete expressions are presented in Appendix B. We interpret this 



as an instability of the deconfined phase in the original theory, since the instantons are 
always in the plasma phase, and the interaction between two well-separated instantons is 
exponentially screened. 

Previous renormalizat ion group analyses of instanton effects in U(l) spin liquids by Na- 
gaosa |10||, Ioffe and Larkin [26(, and Wen accounted only for the flow of the instanton 
fugacity in Eq. (fT9~j). Taken on its own, this equation would indeed suggest that there is a 



regime of parameters (87r ya^a^ — (l/2£) 2 < 1/2 for Eq. (fT9])) where the fugacity y flows to 
zero, and the instanton effects are negligible: this was the conclusion of these earlier works 
However, this conclusion does not hold after accounting for the flows in Eqs. 

(gZHHD - 

IV. CONCLUSIONS 

The primary conclusion of this paper is that a two (spatial) dimensional system of 
fermionic spinons residing on an incipient Fermi surface, and interacting with a compact 
U(l) gauge force, is generically in a confining phase at zero temperature. Consequently, 
there are no free spinon-like excitations. When combined with earlier results on the confine- 
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31 



32, 



our results rule out the 



ment of spinons with a Dirac fermion spectrum 
existence of a spin liquid phase in two spatial dimensions with deconfined spinons interact- 
ing with a U(l) gauge force. In particular, it appears that none of the algebraic, U(l), and 
SU(2) symmetric spin liquids in Wen's classification 1(J are stable phases of matter. 

Our analysis also points to the existence of possible (multi)-critical points at which the 
confinement picture does not directly hold. These correspond to points where accidental 
cancellations cause the renormalized values of a& and/or a w to vanish. These are likely 
quantum critical points between different confining phases, and not conventional Coulomb 
states. Similar points appeared in analyses of quantum dimer models j3, Q • 

Finally, we note that the arguments of this paper do not extend to three spatial dimen- 



sions. Here, a compact U(_ 
spin liquids are possible [4j 



) gauge field can be in a Coulomb phase, and interesting U(l) 
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APPENDIX A: SCREENING OF INSTANTON INTERACTIONS 



The analysis in the body of the paper, and in the Ref. 33 has been carried out using 
the dual sine-Gordon representation. While this is technically a useful way to proceed, it 
does have the disadvantage of obscuring some of the underlying physics. To remedy this, 
we present our arguments here using the direct partition function for the instantons. For 
simplicity, we will restrict our discussion here to the case of Dirac fermion spinons, where the 
'relativistic' nature of the theory does simplify the formalism by making the physics isotropic 



in spacetime. However, similar arguments can be advanced a 



The computations be 



ow expand on the remarks in Refs. 



29 



so for t 



re non-relativistic case. 



30 



31 



32l and were implicit! 



ciny 

y 



carried out in Ref. [n. They also connect with some of the arguments presented in Ref 
as argued below. 

With spinons with a Dirac fermion spectrum, two instantons have an action which de 



271 ] . Let us consider a model 



pends logarithmically on the separation between them [la, |2 
of integer 'charges' rrii (instantons) at spacetime positions Ri, with the interaction term in 
the action J2i<j m i m jV(Ri ~ Rj), where 



V(R) = -K\n{\R\/a) 



(Al) 



with K a coupling constant. We consider such an interaction in D spacetime dimensions. 
We are interested here is D = 3, but the general D analysis will also allow us to connect 
with classical results on the Kosterlitz-Thouless transition. 

Following the original analysis of Kosterlitz j^, we compute the effective interaction 
V e g(R) between two well separated test charges: a +1 charge at Ri and a —1 charge at R2. 
To lowest order in the instanton fugacity, y, this interaction is renormalized by instanton 
dipoles. In the spirit of the renormalization group, we consider only the renormalization by 
dipoles consisting of opposite charges with spacetime separation between a and a + da. So 



— * — * 

the dipole consists of a +1 charge at R + a/2 and a —1 charge at R — a/2. As in Ref. |48, 
to leading order in y the expression for V e s is 

exp (yiRi - R 2 



exp (y eS (Ri - R 2 



1 + 



ra-\-da f —> / _» _» _» 

y 2 J d D a J d D R exp (-V (Hi - i? - a/2) + - i? + a/2) 



+V(H 2 - H - a/2) - V{R 2 -R + a/2)) + £>(y 



where y is the charge fugacity, and the normalization 

ra+da 



d D a I d u R + 0(y 4 



D 



(A2) 



(A3) 



The integral in Z is proportional to the volume of the system: however, this is standard in 
the low density expansion of any system, and the volume dependence cancels in all physical 
quantities order-by-order in y; the logarithm of the partition function is extensive, and this 
dependence on the system volume is merely a signal of that. In the above expression for 
V e ff it is assumed that any two charges are never less than a distance a apart, and this cuts 
off any incipient ultraviolet divergences. Now expand the argument of the exponential in 
powers of a. This yields 

exp (V(Ri - R 2 ) 



exp [V eS (R 1 - R 2 



1 + 



y 2 f a+da d D d I d D Rexp -Kd- 1 ' h " 



Rn — R 



0(d 3 ) + 0(y 4 



(i?i - R) 2 (R 2 -R) 2 
The integral over a can be performed, and to leading order in y and a we obtain 



(A4) 



Vcs(Ri — R 2 ) 
= -Kln(\R 1 -R 2 \/a)+y 

= -Kln(\Rx-R 2 \/a)+y 



2 S D K 2 a D+1 da 



2D 



d D R 



Ri — R 



Rn — R 



2 S D K 2 a D+l da 



2D 



d D R- 



{Ri-R) 2 (R 2 -R) 2 

(Ri - R2) 2 
(Ri - R) 2 {R 2 - R) 2 



(A5) 



(A6) 



where is the surface area of a sphere in D dimensions. Note that the volume dependence 
has cancelled as expected. The above is, of course, an adaptation of the computation 
performed by Kosterlitz j^sj ]. 

Now, notice that there is a crucial distinction in the integral in (|A6|) between D = 2 
and D — 3. In D — 2, the integral reproduces the functional form of the bare logarithmic 
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interaction, and so the renormalization of the potential can indeed be accounted for by a 
dielectric constant: 

- - S 2 K 2 o?da -> 
V eS (R 1 -R 2 ) = -K\n(\R 1 -R 2 \/a)+y 2 ^-^ [H\Ri - R 2 \ /a) + constant] , D = 2, 

(A7) 

where the ultraviolet cutoff radius around each charge is used to make the integral finite. 
This leads to the usual Kosterlitz recursion relations |4£| in D = 2. 
Turning, finally, to D = 3, the integral over R in (jA6|) now yields 



- ^ 2 TT 2 K 2 n 4 dn -> 

V eS (R 1 -R 2 ) = -Kln(\R 1 -R 2 \/a)+y 2 ^^ \Ri - R 2 \ , D = 3. (A8) 

There was now no need to apply any ultraviolet cutoff. Note that the bare logarithmic 
interaction is not reproduced under renormalization-this is the central point. Instead, we 
have generated a linear interaction, which has a much stronger dependence on the separation 
between the charges. It is also crucial to note the sign of this linear interaction: it is such as to 
weaken the interaction between oppositely charged particles, and at large enough R\ — R 2 
the effective interaction between these opposite charges actually becomes repulsive. Of 
course, we cannot trust the renormalization group beyond these scales. Nevertheless, the 
lesson is quite clear: the effect of the dipole is to substantially weaken the interaction between 
the charges, suggesting a ubiquitous instability to a plasma phase. 

Continuing this line of argument, we can also look at the screening of the generated linear 
interaction. Using exactly the same derivation as above, the integral in the second term on 
the right hand side of (|A5|) will be replaced by 

d 3 n ( 4^4 - 4^4 V (A9) 

\ \Ri — R\ | -R2 — R\ J 

This integral is now divergent at large R, which reconfirms our conclusion that the screening 
effect of the largest scales R is extremely disruptive. 

TheEq. (jSHJ) may also be understood as the expansion of the screened interaction found 
in Ref. |3J2 in small fugacity. Considering the electrostatics of logarithmically interacting 
charges bound in dipole pairs, the effective interaction due to an external charge, in the 
Fourier space, becomes, 

Kfr(g) = | ,3 f 2 , (Aio) 



if 



in D = 3, where x ~ y 2 is the finite polarizability of the medium 33]. At a large distance R 
this translates into V e ff(R) ~ l/(47ixR). Expanding in small \i on t ne other hand, yields 

K Ky 

V cS (q) = -- 3 -^ + 0(x 2 ). (All) 
\Q\ q 

The first term in the expansion then represents the bare logarithmic interaction in Eq. (|A8|) . 
while the second one may be identified with the next, linear correction. 



APPENDIX B: COMPUTATION FOR RENORMALIZATION ANALYSIS 

Here we discuss the smooth cutoff procedure we devised to compute G > (x, r) in Eq. (fTtjj). 
and to arrive at the flow equations (jl7H19|) . With a u = in the unrenormalized theory, 

G>(x,r) = — ^ / —J (kx) / — ; 2 , v Bl 

First, we write 

/ A ^ J (A;x) = r ^ J (fcr) - /°° ^J (kx) (B2) 
fc | — JA/b k Ja k 

and then introduce a smooth cutoff |33J] as 

A T Jo(te) ^7o I J ° (b fe- (B3) 

This leads to 

1 f A dk T n xAln(b) dK (z) . 

with i^o(^) as the Bessel function. Similarly, for Q/e 2 £ C 1 we may neglect the frequency 
dependence and approximate the second integral in Eq. (|B1|) as 

° r d ^r e -^, m ? (B5) 

-a 2tt 7-oo 2tt ' v ; 

where we have introduced the Gaussian instead of the sharp cutoff in the frequency space. 

This gives 

n d_^^z y n (TO)a/ir 

2tt 1 + |w|/(e 2 £) vr { ' 

Inserting the last two expressions into Eqs. (fT3*|) and (fPfj) yields the ~ y 2 terms in the 

recursion relations Eqs. (fTTj) and (fTHj). 

Assuming > 0, the frequency cutoff f2 may be taken to infinity in the computation of 
G>(0, 0) appearing in the Eq. (|T5j). This leads to 



dln(fo) 



2 / 1 

1 arctan 
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y + 0(f), (B7) 



for D > 0, where 



D = ia k aj 2 - 1. 



(B8) 



For D < 0, one finds 



dy 



d\u(b) 



2 - 



4^2 mi 1/2 



hi 



1 + 


D 


1/2 


1 - 





(B9) 



For completeness, let us also mention the case a w (b) = 0, which would arise if one would 
neglect ~ y 2 terms in the recursion relations. One then needs to keep a finite cutoff Q, to 
find 



dy 



d\n(b) 



2 - 



4vr 2 



In 1 + 



y + o(f). 



(bio) 



Because to the lowest order in y neither £ nor dk scale, the last equation would suggest that 
fugacity is irrelevant for small enough d k , for example. This conclusion, however, inevitably 
breaks down to the next order in y, as discussed in the body of the paper. 

Note also the similarity between the Eq. (|B10)1 and Eq. (11) in Ref.llO, upon identification 
of i with 7, and a k with g 2 . The main difference is that in Nagaosa's renormalization scheme 
the coupling g 2 (analogous to our coefficient a^) appears irrelevant by naive power counting, 
which then suggest that fugacity should become irrelevant as well. This conclusion, however, 
changes to the next order in fugacity, in which the fugacity always becomes relevant even 
within Nagaosa's renormalization scheme, and much like in our calculation. 
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